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NONLINEAR  THEORY  OF  THE  E  X  B  INSTABILITY  WITH  AN 
INHOMOGENEOUS  ELECTRIC  FIELD 


I.  INTRODUCTION 

The  E_  x  _B  instability,  also  known  as  the  gradient-drift  instability, 

has  been  invoked  to  explain  both  natural  and  artificially  induced  plasma 

density  structure  and  irregularities  in  the  terrestrial  ionosphere.  This 

instability  can  be  excited  in  a  low  pressure,  weakly  ionized,  magnetized 

plasma  that  contains  ait  ambient  electric  field  orthogonal  to  both  a 

magnetic  field  and  a  density  gradient.  The  basic  physical  mechanism  for 

the  E_  x  _B_  instability  [Linson  and  Workman,  1970;  Perkins  et  al. ,  1973]  is 

analogous  to  that  describing  the  classical  Rayleigh-Taylor  instability  in 

which  a  heavy  fluid  is  supported  by  a  lighter  fluid.  Originally  applied  by 

Simon  [1963]  and  Hoh  [1963]  to  laboratory  gas  discharges,  this  instability 

has  been  applied  to  ionospheric  plasma  cloud  structuring  [Haerendel  et  al. , 

•• 

1967;  Linson  and  Workman,  1970;  Volk  and  Haerendel,  1971;  Perkins  et  al.  , 
1973;  Zabusky  et  al. ,  1973;  Scannapieco  et  al. ,  1976;  Ossakow  et  al. ,  1978; 
Chaturvedi  and  Ossakow,  1979;  Keskinen  et  al. ,  1980;  McDonald  et  al.,  1981] 
and  to  the  stability  and  transport  of  large  scale  convecting  auroral 
ionospheric  plasma  enhancements  [Keskinen  and  Ossakow,  1982,  1983;  Vickrey 
et  al. ,  1980].  However,  these  studies  have  addressed  only  the  E  x  B 

instability  driven  by  initial  ambient  homogeneous  electric  field.  Perkins 
and  Doles  [1975]  have  shown,  using  linear  theory,  that  sheared  velocity 
flow  (resulting  from  an  initially  self-consistent  inhomogeneous  electric 
field  parallel  to  the  density  gradient)  can  stabilize  the  _E_  x  B_  instability 
in  the  collisional  regime.  Huba  et  al.  [1983]  verified  numerically  the 
results  of  Perkins  and  Doles  [1975]  and  extended  the  linear  theory  of  the  _E_ 
Manuscript  approved  October  3,  1983. 
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x  B  instability  with  an  inhomogeneous  electric  field  to  the  collisionless 
regime.  In  order  to  properly  apply  the  E_  x  3_  instability  to  the  auroral 
and  polar  ionosphere  and  magnetosphere  one  must  consider  the  linear  and 
nonlinear  evolution  of  the  E  x  B  instability  with  an  inhomogeneous  electric 
field  since  it  is  well  known  that  high  latitude  magnetospheric  and 
ionospheric  electric  fields  are  usually  inhomogeneous  [see,  for  example, 
Fairfield,  1977  and  references  therein]. 

The  purpose  of  this  paper  is  to  study  the  nonlinear  evolution  of  the 
E_  x  B_  instability  with  an  inhomogeneous  electric  field.  The  organization 
of  the  paper  is  as  follows.  In  Section  2  we  give  the  general  equations 
describing  the  E  x  B  instability  with  an  inhomogeneous  electric  field.  In 
Section  3  we  study  the  nonlinear  evolution  of  the  _S_  x  _3_  instability  with 
inhomogeneous  electric  field  by  numerically  solving  the  fundamental  fluid 


& 


+1* 


equations.  In  Section  4  we  give  a  nonlinear  analytic  theory  of  the  results 
in  Section  3.  Finally,  in  Section  5  we  discuss  and  summarize  our  findings. 


2.  MODEL  EQUATIONS  AND  LINEAR  THEORY 

The  basic  equations  describing  the  evolution  of  the  _E_  x  _B_  instability 
in  an  inhomogeneous  electric  field  are: 


9N 


— -  +  V.fN  V  )  -  0 
3t  ra 


-  —  fE  +  c_1V  x  B]  -  0 

m  — e  — ' 


s-  (E  +  *  B)  -  vln  -  0 


V.J  -  7*n(v_1  -  Vj  -  0 


(1) 


(2) 


(3) 


(4) 


where  a  denotes  species  (ct  ■  e,i)  and  other  symbols  retain  their 
conventional  meaning. 


The  equilibrium  configuration  used  in  the  analysis  is  shown  in  Fig. 
1.  The  ambient  magnetic  and  electric  fields  are  in  the  z  direction  and  the 
x,y  plane,  respectively,  where  B  -  and  E  =  Eqx(x)  x  +  E  y_.  The 
electric  field  in  the  y  direction  is  constant,  while  the  electric  field  in 
the  x  direction  is  allowed  to  be  a  function  of  x.  This  gives  rise  to  an 
inhomogeneous  velocity  flow  in  the  y  direction,  i.e.,  7Qy(x)  M 
-cEQX(x)/B.  The  density  is  taken  to  be  inhomogeneous  in  the  x  direction 
(nQ  *  nQ(x))  and  temperature  effects  are  ignored. 

The  basic  assumptions  used  in  the  analysis  are  as  follows.  We  assume 
that  the  perturbed  quantities  vary  as  of  -  5f(x)  exp[i(k^y  -  jt'l,  where  k^ 
is  the  wave  number  along  y  direction  and  a  ■  jJr  +  iy,  implying  growth  for 
y  >  0.  The  ordering  in  the  frequencies  is  such  that  u>  «  ^  and  << 

(F  region  approximation),  where  vin  is  the  ion-neutral  collision  fre-quencv 
and  is  the  ion  gyrof requency .  We  ignore  finite  gyroradius  effects  by 
limiting  the  wavelength  domain  to  krLi  <<  1,  where  rLi  is  the  mean  ion 
Larmor  radius.  We  neglect  perturbations  along  the  magnetic  field 
fk  »  0)  so  that  only  the  two-dimensional  structure  in  the  x,y  plane  is 
obtained. 

Solving  Eq.  (2)  and  (3)  for  and  V^,  substituting  into  (1)  (for 
electrons)  and  (4)  we  find  [Perkins  and  Doles,  1975] 


Perking  and  Doles  (1975)  expand  Eq.  (8)  about  x  •  xQ  where  x0  is  the 
position  of  maximum  n'/nQ  by  taking  (f'  =  df/dx) 

n  '/n  -  ri-fx-x  U/D2VL  (9) 

o  o  >■  '  o'  ■  n 
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Assuming  k^L^  >>  1  and  k2D2  <<  1,  and  by  making  several  variable  changes, 
they  solve  Eq.  (3)  analytically.  The  important  conclusion  of  their  results 
is  that  the  E_  x  3_ instability  is  stabilized,  in  the  linear  regime,  when 


E  (x  ) 
x  o  v  2 

E  >  k~D 

y  y 


(10) 


Thus,  the  influence  of  velocity  shear,  i.e.,  an  inhomogeneous  Eox,  is  to 
preferentially  stabilize  the  short  wavelength  modes,  i.e.,  those  with  kvD 

j 

>>  1.  Huba  et  al.  [1983]  verified  Eq.  (10)  by  solving  numerically  the 
fully  nonlocal  equation  (8). 

To  recover  the  local  theory  limit  we  let  ->•  ik  and  assume 

3  X  x 

k^L^  »  1  and  k^L^  »  \  where  •  (n'/n^)  ^  is  the  scale  length  of  the 
density  inhomogeneity  evaluated  at  x  ■  xQ.  For  simplicity  we  also  take 
_E  *  Eqx  x  +  EQy  jy  *  constant.  Following  standard  techniques,  we  find  for 
the  frequency  and  growth  rate 
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k  V 
y  oy 
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(11) 


(12) 


the  usual  _E_  x  _B_  instability  growth  rate  [Linson  and  Workman,  1970]  in  the 


collisional  limit. 


3. 


NONLINEAR  EVOLUTION 


Since  the  E_  x  _3_  instability  becomes  highly  nonlinear  and  analytically 
intractable  we  will  study  the  nonlinear  evolution  of  the  E_  x  3_  instability 
in  an  inhomogeneous  electric  field  by  numerically  solving  the  fundamental 
equations  (5")  and  (6).  We  choose  parameters  typical  of  naturally  occurring 
high  latitude  convecting  ionospheric  plasma  enhancements  [Vickrey  et  al.  , 

1980]  and  artificially  produced  ionospheric  plasma  clouds  [McDonald  et  al.  , 

1981] .  Equations  (5)  and  (6)  were  solved  on  a  numerical  grid  consisting  of 
258  grid  points  in  the  x-direction  and  102  grid  points  in  the  y-direction 
with  constand  grid  spacing  of  0.3  km  in  the  x-direction  and  0.2  km  in  the 
y-direction.  As  a  result,  the  simulation  plane  has  an  x,y  extent  of  30  and 
20  km,  respectively.  The  plasma  density  N  in  equation  (5)  was  advanced  in 
time  using  a  multi-dimensional  flux-corrected  variable  timestep  leapfrog- 
trapezoid  scheme  [Zalesak,  1979]  which  is  second  order  in  time  and  fourth 
order  in  space.  At  each  timestep  the  self-consistent  electrostatic 
potential  5<j)  of  the  plasma  enhancement  in  Eq.  (6)  was  deter-mined  using  a 
Chebychev  iterative  method  [McDonald,  1980]  with  a  convergence  criterion  of 
10~4.  Periodic  boundary  conditions  were  imposed  in  the  y-direction  with 
Neumann  boundary  conditions  (3/3x  *  0)  in  the  x-direction.  A  slab 
approximation  is  used  to  model  the  zeroth  order  plasma  density  with  profile 
given  by  nQ(x)  *  NQ  {1  +  4.5  [tanh(x-x^)/L^  +  tanh  (X2~x)/L>j]’  ,'l  +  c(x,y)' 
with  x^  *  10  km,  X2  *  35  km,  and  L^  =  10  km.  This  gives  a  maximum  plasma 
density  to  background  ratio  of  approximately  10.  The  initial  perturbation 
in  taken  to  be  completely  random  with  root-mean-square  amplitude  of  10-4. 
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The  ambient  electric  field  is  chosen  to  be 


E  (x) 


E  (x) 

ox 


X  + 


(13) 


where 


E  (x  ■  -  ®)  »  E  sin  9  x  +  E  cos  9  y  (14) 

~0  O  —  0  — 

so  that  9  »  tan“^(E0X/E0y)  at  x  3  The  influence  of  the  x  component  of 

the  electric  field  is  then  studied  by  varying  9,  the  angle  between  E  and  _v 
at  x  ■  -  <».  The  form  of  EQX(x)  considered  in  the  analysis  is: 

Eqx(x)  -  Eq  sin  9(NQ/no(x))  t  constant  (15) 

We  comment  that  Eq.  (15)  is  an  equilibrium  solution  which  satisfies 
V*.J  ■  0  i.e.,  Eq.  (4)  or  (6). 

We  consider  two  models  with  different  initial  electric  field 
configurations  to  illustrate  the  effect  of  an  inhomogeneous  electric 
field.  Model  1  has  EQy  *  10  mV/m,  EQX  ■  0  (no  velocity  shear)  while  Model 
2  takes  Eox(-  «)  «  8.6  mV/m  and  EQy  *  5  mV/m  giving  9  =  60°.  These 
electric  field  magnitudes  are  chosen  to  be  typical  of  the  high  latitude 
diffuse  auroral  F-region  ionosphere  [Vickrey  et  al. ,  1980].  Figure  2a-2d 
shows  the  evolutions  of  the  £  x  B  instability  using  Model  1  (no  velocity 
shear).  Figure  2a  shows  the  initial  configuration  which  includes  the  small 
random  perturbation.  Figure  2b  illustrates  the  linear  regime  at  t  =  250 
sec  (yt  *  5)  and  shows  unstable  growth  on  the  trailing  side  of  the  plasma 
enhancement  as  predicted  by  the  linear  result  given  by  Eq.  (12).  One  can 


note  the  depletion  jetting  to  the  front  side  of  the  enhancement  in  analogy 
to  the  initial  evolution  of  the  _E_  x  3_  gradient  drift  instability  in 
artificial  ionospheric  plasma  clouds  [Zabusky  et  al.  ,  1973  ;  Scannapiecoet 
al. ,  1976].  Figure  2c  gives  the  structure  of  the  plasma  enhancement  at  t  = 
500  sec  and  shows  steepened  fingers  which  are  beginning  to  elongate. 
Finally  Figure  2d  displays  the  plasma  enhancement  at  t  =  650  sec  in  the 
fully  nonlinear  regime.  The  trailing  edges  of  the  principal  fingers 
(striations)  have  steepened,  become  quasi-one  dimensional  and  bifurcated. 
The  length  scales  on  Figure  2a-d  are  distorted  with  the  depletions  longer 
and  narrower  than  is  depicted. 

Figure  3a-d  give  sample  one-dimensional  spatial  power  spectra  at  t  ■ 
0,  250,  650  sec  Model  1.  These  power  spectra  are  defined  as  follows 

?<v  » ;  ?(kx,ky) 

p(ky)  =  ;  dkx  p(kx,ky) 

_  2 

where  P(k  ,k  )  =  (L  L  )  r£n(k  ,k  )/N  1  is  the  spectral  density,  on  ■  N-N 
x  y  x  y  -  x  y  o  1  o 

with  NQ  the  peak  plasma  enhancement  density,  and  L^Ly  is  the  area  of  the 
numerical  simulation  plane.  Figure  3a  illustrates  the  power  spectrum 
(P(ky)  of  the  random  perturbation  used  to  initialize  the  Z_  x  3_ 
instability.  Figure  3b  shows  the  power  spectrum  P(ky)  in  the  linear  stage 
of  the  instability  and  compares  favorably  with  the  growth  rate  vs.  k,,  from 


local  theory  [Huba  et  al. ,  1983].  Figure  3c-3d  gives  the  power  spectra  in 
the  x-  and  y-directions ,  respectively,  in  the  nonlinear  regime  at  t  =  650 


sec.  For  both  cases  these  power  spectra  are  well-fitted  with  an  inverse 
power  law  with  spectral  index  nx  =  2  for  2-r/k^  between  approximately  30  and 
1  km  and  n„  a  2-2.5  for  2-rr/k  between  20  and  1  km. 

y  y 

Figures  4a-4d  illustrate  the  evolution  of  the  _E_  x  _B_  instability  with 
an  inhomogeneous  electric  field  using  Model  2.  Figure  4a  gives  the  initial 
configuration  which  is  identical  to  Fig.  2a.  Figure  4b  shows  the 
isodensity  contours  of  the  plasma  enhancement  in  the  linear  regime  at  t  = 
700  sec  (yt  a  5).  Figure  4c  gives  the  evolution  at  t  =  900  sec  and  shows  a 
distinct  bending  of  the  striations.  Figure  4d  details  the  E  x  3 
instability  in  the  nonlinear  regime  at  t  =  1250  sec  where  one  notes  that 
the  fingers  (striations)  are  no  longer  primarily  aligned  with  the  flow. 
Similar  morphologies  are  also  observed  for  other  velocity  shears 
(9  -  30°). 

Figures  5a-5c  give  sample  power  spectra  for  Model  2  in  the  linear  and 

nonlinear  regime  at  t  ■  425  and  1250  sec,  respectively.  Figure  5a 

illustrates  the  linear  regime,  shows  the  suppression  of  the  shorter 

wavelength  fluctuations  in  agreement  with  linear  theory,  eq.  (10),  [ Perkins 

and  Doles,  1975;  Huba  et  al. ,  1983],  and  indicates  a  preferred  scale 

size.  Figure  5b  gives  a  sample  power  spectrum  in  the  y-direction,  P(ky), 

-n 

and  can  be  described  by  a  power  law  P(k^)  -  k^  ,  n  =  2-2.5  in 

approximate  agreement  with  the  no  shear  case  Model  1.  The  power  spectrum 

in  the  x-direction  P(kx)  in  the  nonlinear  regime  is  given  in  Fig.  5c  and 

-n  _ 

can  also  be  described  by  a  power  law  P(kx)  ~  k  ,  n^  =  2.  Similar  power 
laws  and  spectral  indices  are  also  observed  for  the  case  where  9  =  30°. 

Figure  6  shows  the  evolution  of  the  mean  density  profile  along  the  x- 
axis  (averaged  over  the  y-direction)  at  t  =  0,  700,  950  sec  during  the 
linear  unstable  stage  of  Model  2.  Similar  result  are  found  for  the  case 


r-  -  .  r.  ». 


3  ■  30°.  We  show  in  the  next  section  that  this  relaxation  is  responsible, 
in  part,  for  the  stabilization  of  the  fastest  growing  linear  nodes. 

4.  QUASILINE AR  THEORY 

We  present  arguments  for  a  quasilinear  stabilization  mechanism  for  the 
E_  x  B_  instability  with  an  inhomogeneous  electric  field.  We  show  that  the 
mean  density  gradient  driving  the  instability  is  modified  by  an  unstable 
spectrum  of  waves  of  finite  amplitude.  By  dividing  N  =  nQ  +  on,  + 

+  6V<_ ,  <5n>  *  <SVe>  *  int0  mean  and  oscillating  parts  (with  <  >  denoting 
space  and  time  average),  the  electron  continuity  equation  (1)  can  be 
written 


5no  =  3<N> 

at  =  at 


-  7*<6n  oV  > 


(16) 


+  n  V •  <5 V  +  5n7*V  +  V  .75n  +  5V  .7n 

at  o  — e  — eo  — eo  — e  o 


7*<5n  5V  >  -  7*Sn  5V 


(17) 


with  <5n  5V  >  »  (L  )_1  f  dy  5n  <$V  -  (L  )_1  f 

— e  y  '  —  — e  y 


5V 


dk 

77^3  5Ve,-k 

(2tt)  -y  -y 


-  -jj-  75*  x  2  and  Ly  is  the  length  of  the  system  in  the  y-direction. 


To  find  3nQ/at  in  Eq.  (16)  to  lowest  (quadratic)  order  in  5®  one  needs  to 
compute  5n  to  linear  order  in  5<|i  using  Eq.  (17)  which  gives: 


Snk  (x) 


-  ^  k_  ^  Fu-k  V_.(x)1_154)  (x) 


c 

3 


y  3x 


y  ey' 


^k 

-y 


(18) 
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<5n  5V  >  -  (L  )  k  x  z  r 

— e  y  3  -y 


- —  5*-k.  5nk 

(  2-r  )  3  -y  -y 


with  5n^  given  by  (18).  Inserting  into  (16)  we  find 


o  3  .  3 

— — •  ■  -  — —  D(x)  —  n 
3t  3x3xo 


2  ky\  (t) 

with  D(x)  -  Re  [i  /  dky  ^_k  vy-  -fx-; 


=  —  f  dk 


!k  (O 


B2  '  7  fu.  “k V_(x)12+yS 


-  r  eo 


where  a  ■  ojr  +  iy  has  been  used.  The  spectral  energy  density 
Ik  (t)  2  ( t )  | 2  changes  with  time  according  to 


3lk  (t) 


2\(t)Ik  (t) 

-  -y 


Equation  (19)— (21)  provide  a  description  of  the  quasilinear  evolution  of 

the _E_ x  B_ instability  with  an  inhomogeneous  electric  field. 

Assuming  weakly  growing  modes,  i.e.,  «  Im  ~  kV  (x)|2  the 

ic  r  <20 

diffusion  can  be  approximated  by 


D<*>  l  dk  Y^k^  <t). 

-y  - y 


From  the  work  of  Huba  et  al.  [1983]  and  the  previous  numerical  simulation 


11 


results  we  observe  that,  although  y^  maximizes  for  preferred  scale  size 

- y 

given  by  k^L  -  0(1),  the  distribution  of  y  vs.  ky  is  broad.  As  a  result, 

we  take  y,  to  be  of  the  fora 
k 

-y 

\  (t)  "  Y0(C)exp^(kv  ‘  koi2/kw: 

-y 

In  addition,  we  assume 


\  (t)  -  Io(t)exp[-(k  -  \)2/'*-2} 

-y 


,-i 


where  yQ  =  £Vq/L,  kQ  =  qL  ,  and  is  the  mean  width  of  the  distribution 
with  kw  >  kQ.  Here,  s,q  are  constants  of  order  unity  [Huba  et  al.  ,  1983. 


Inserting  these  expressions  for  y ,  and  1^  into  eq.  (22)  we  take 

iy  -y 

k 


w 


D(x)  =  ll  y-1I  f  °  l  dk  k2 
B2  o  o  '  k  y  y 

k - - 

o  2 


c2  L 

= 

s  o 


k3 

w 

[  — — 
o  12 


(23) 


For  approximate  nonlinear  saturation  of  the  fastest  growing  mode  with 

wavenumber  k  we  have 
max 


Y 


yt  -  Dk2 
•L  max 


(24) 


Using  eq.  (23)  in  (24)  we  find 


where  <j>  «  (BV  L/c).  From  Huba  et  al.  [19831,  :  3  0.2,  k  L  =  1, 

Yo  o  -  L  max  * 

k  L  »  5,  which  gives  5<s/i»  ~  0.11  which  is  not  inconsistent  with  the 

w  o 

previous  numerical  simulation  results  in  Sec.  3.  After  the  initial  mean 
density  gradient  has  relaxed  to  a  certain  degree,  small  scale  modes  (kL  >> 
1)  can  become  unstable  leading  to  two-dimensional  wave  coupling  processes 
which  will  determine  the  wave  spectrum. 

5.  DISCUSSION  AND  SUMMARY 

We  have  studied,  using  analytical  and  numerical  techniques,  the 
nonlinear  evolution  of  the  E_  x  B_ instability  with  an  inhomogeneous  electric 
field.  The  principal  results  of  this  study  are  as  follows: 

1.  The  basic  morphology  and  power  spectra  of  the  E_  x  3_  instability 
with  an  inhomogeneous  electric  field  in  the  nonlinear  regime  is 
similar  to  homogeneous  (no  shear)  case. 

2.  Nonlinear  effects  destabilize  the  linearly  stable  modes  associated 
with  the  E_ x  B_ instability  with  inhomogeneous  electric  field. 

3.  Quasilinear  mechanisms,  in  which  the  initial  density  gradient 
driving  the  inhomogeneous  E_  x  B_  instability  is  modified  by  a 
finite  amplitude  wave  spectrum,  can  contribute  to  nonlinear 
stabilization. 

.  Recently,  large  scale  equatorward  convecting  plasma  enhancements  in 
the  diffuse  auroral  F-region  ionosphere  have  been  identified  and  studied 
[Vickrey  et  al. ,  1980]  using  both  radar  and  satellite  measurements. 
Observed  in  regions  of  diffuse  auroral  particle  precipitation  and 
associated  field  aligned  currents,  these  enhancements  have  overall 
latitudinal  dimensions  of  a  few  hundred  kilometers,  contain  relatively 


steep  poleward  and  equatorward  edges,  and  have  been  shown  to  be 
approximately  field-aligned  resembling  vertical  slabs  of  ionization.  Their 
occurrence,  which  is  maximized  in  the  evening-midnight  sector,  is 
apparently  not  strongly  related  to  magnetic  activity  nor  to  Z-region 
processes.  The  presence  of  plasma  density  irregularities  associated  with 
these  enhancements  has  been  verified  using  satellite  scintillation  studies 
[Fremouw  et  al.  ,  1977;  kino  et  al.  ,  1978;  Vickrey  et  al.  ,  19801.  The 


scintillation  data  have  indicated  that  the  electron  density  irregularities 
are  structured  like  L-shell  aligned  sheets  for  irregularity  scale  sizes  of 
approximately  1  km  [Fremouw  et  al.,  1977;  kino  et  al.  ,  1980].  Moreover, 
the  source  region  of  these  scintillation  causing  irregularities  has  been 
demonstrated  to  be  latitude  limited  [kino  and  Owen,  1980]  and  contained  in 
a  vertical  slab  of  F  region  plasma.  Keskinen  and  Ossakow  [1982,  1983] 
showed  that  these  convecting  plasma  enhancements  can  be  destabilized  by  the 
E  x_B  instability.  Their  results,  using  a  homogeneous  electric  field,  were 
not  inconsistent  with  available  experimental  observations  [kino  et  al. , 
1978;  Vickrey  et  al. ,  1980;  Tsunoda  and  Vickrey,  1982].  The  results  of  the 
present  study,  with  an  inhomogeneous  electric  field,  may  help  explain  the 
geometric  (L-shell  alignment)  nature  of  these  structures  if  one  assumes  the 
east-west  (north-south)  direction  corresponds  to  the  y(x)-axis.  In  this 
case,  the  fingerlike  structures  in  the  nonlinear  regime  of  Model  2  would  be 
approximately  east-west  or  L-shell  aligned. 

These  results  are  also  applicable  to  the  development  of  the  _E_  x  _3 
instability  in  artificially  injected  ionospheric  plasma  (barium)  cloud 
releases.  3y  including  a  self-consistent  inhomogeneous  electric  field,  we 


note  that  the  time  scales  for  striation  formation  and  jetting  is  increased 


compared  to  the  homogeneous  case.  This  effect  will  lead  to  an  increase  in 
striation  onset  time  [McDonald  et  al.,  1981].  However,  in  the  nonlinear 
regime,  our  studies  with  an  inhomogeneous  electric  field  show  similar 
morphologies  and  power  spectra  with  respect  to  the  homogeneous  case. 
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(C)  (d) 

Fig.  2  —  Real  space  isodensity  contour  plots  of  n(x,y)/N0  for  model  1  at  (a)  t  *  0  sec, 
(b)  t  *  250  sec,  (c)  t  ■  500  sec,  (d)  t  =  650  sec.  The  x-axis  is  compressed  by  a  factor  of  2.58. 
The  distance  between  tic  marks  in  the  x-direction  (y-direction)  is  5  km  (12.8  km).  Eight 
contours  are  plotted  in  equal  increments  of  1.25  beginning  at  1.25.  The  observer  is  looking 
downward  along  the  magnetic  field  lines. 
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Fig.  3  —  Sample  one-dimensional  power  spectra  for  model  1  at  (a)  t  *  0  sec,  P(ky)  (b)  t  = 
250  sec,  P(ky)  (c)  t  *  650  sec,  P(kx)  (d)  t  =  650  sec,  P(ky).  Here  kpx  =  2^/80  km*l  and 
kpy  *  2jt/20  km-1,  the  solid  curve  is  a  least  squares  fit  to  modes  2-80  in  the  x -direction  and 
to  modes  2-30  in  the  y -direction.  The  units  of  P(kx ),  P(kv )  are  km. 
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Fig.  5  —  Sample  one-dimensional  power  spectra  for  model  2  at  (a)  t  =  425 
(b)  t  =  1250  sec,  (c)  t  =  1250  sec  in  same  format  as  Figure  3. 
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ARLINGTON,  VA  22217 
01CY  ATTN  CODE  465 
01CY  ATTN  CODE  461 
01CY  ATTN  CODE  402 
01CY  ATTN  CODE  420 
01CY  ATTN  CODE  421 


AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  NM  87117 
01CY  ATTN  SUL 

01CY  ATTN  CA  ARTHUR  H.  GUENTHER 
01CY  ATTN  NTYCE  1LT.  G.  KRAJEI 

AFTAC 

PATRICK  AFB,  FL  32925 
01CY  ATTN  TF/MAJ  WILEY 
01CY  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  AAD  WADE  HUNT 
OICY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  &  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
OICY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 
OICY  ATTN  J.  DEAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
OICY  ATTN  YSEA 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVIS’.ON/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  APB,  MA  01731 

01CY  ATTN  DCKC  MAJ  J.C.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  NICD  LIBRARY 

01CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFPISS  AFB,  NY  13441 

01CY  ATTN  DOC  LIBRARY/TSLD 
01CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  ADWATE  MAJ  BRUCE  3 AUER 
01CY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01CY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP  A.  L0RENT2EN 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  309 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

OICY  ATTN  DOC  CON  FOR  J.  WOLCOTT 
OICY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 
OICY  ATTN  DOC  CON  FOR  E.  JONES 
OICY  ATTN  DOC  CON  FOR  J.  MALIK 
OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 
OICY  ATTN  DOC  CON  FOR  J.  ZINN 
OICY  ATTN  DOC  CON  FOR  P.  KEATON 
OICY  ATTN  DOC  CON  FOR  D.  WESTERVELT 
OICY  ATTN  D.  SAPPENFIELD 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUGH 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 
OICY  ATTN  DOC  CON  FOR  T.  COOK 


DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWITZ 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 
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OTHER  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

01CY  (ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

01CY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROMBIE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01CY  ATTN  R.  GRUBB 
01CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01CY  ATTN  I.  GARFUNKEL 
01CY  ATTN  T.  SALMI 
01CY  ATTN  V.  JOSEPHSON 
01CY  ATTN  S.  BOWER 
01CY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

01CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 
01CY  ATTN  L.  SLOAN 
01CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
01CY  ATTN  J.  WORKMAN 
01CY  ATTN  C.  PRETTIE 
01CY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

0 ICY  ATTN  G.  KEISTER 
01CY  ATTN  D.  MURRAY 
01CY  ATTN  G.  HALL 
01CY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
01CY  ATTN  D.B.  COX 
01CY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 
0 ICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
01CY  ATTN  H.  LOGSTON 
01CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wilshire  Blvd. 

Santa  Monica,  Calif  90401 
01CY  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
01CY  ATTN  J.  ROBERTS 
01CY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

01CY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
01CY  ATTN  F.  REIBERT 


GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 
HMES 

COURT  STREET 
SYRACUSE,  MY  13201 
01CY  ATTN  G.  HILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  QIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61320 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  6  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

01CY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  COR? 

P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
0  ICY  ATTN  DASLAC 
OICY  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  T>.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  A VENUS 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CT  ATTN  GEORGE  MROZ 
OICT  ATTN  W.  OLSON 
OICY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
OICY  ATTN  W.F.  CREVIER 
01CY  ATTN  STEVEN  L.  GUTSCHE 
OICY  ATTN  R.  BOGUSCH 
OICY  ATTN  R.  HENDRICK 
OICY  ATTN  RALPH  KILB 
OICY  ATTN  DAVE  SOWLE 
OICY  ATTN  F.  FAJEN 
OICY  ATTN  M.  SCHEIBE 
OICY  ATTN  CONRAD  L.  LONGMIRE 
OICY  ATTN  B.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NEW  MEXICO  37106 
OICY  R.  STSLLINGWERF 
OICY  M.  ALME 
OICY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

OICY  ATTN  JOHN  MORGANSTERN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
OICY  ATTN  W.  HALL 
OICY  ATTN  W,  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LA 3 
313  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LA3 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  6  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 
OICY  ATTN  L.  DeRAND 
OICY  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUCALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
01CY  ATTN  R.  TSUNODA 
01CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  G.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
01CY  ATTN  WALTER  JAYE 
01CY  ATTN  J.  VICKREY 
01CY  ATTN  RAY  L.  LEADABRAND 
01CY  ATTN  G.  CARPENTER 
01CY  ATTN  G.  PRICE 
01CY  ATTN  R.  LIVINGSTON 
01CY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01CY  ATTN  JOHN  ISE,  JR. 

01CY  ATTN  JOEL  GARB AR I JO 


TRW  DEFENSE  4  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01CY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 
OICY  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF 
NOTED) 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

Dr.  P.  MANGE  -  CODE  4101 
Dr.  P.  GOODMAN  -  CODE  4180 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  S WIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  W.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINQZEN 


FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GRESKBELT,  MD  20771 
DR.  K.  MAEDA 
DR.  S.  CURTIS 
DR.  M.  DUBIN 

DR.  N.  MAYNARD  -  CODE  696 
COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 
ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 


IONOSPHERIC  MODELING  DISTRIBUTION  LIS' 
(UNCLASSIFIED  ONLY) 

HE 


INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  E.  FRIEMAN 


STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  3ANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
A3ERDEEN,  MD 

DR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AX  99701 
DR.  L.E.  LEE 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-iO,  MS-664 
LOS  ALAMOS,  NM  3*545 
DR.  M.  PONGRATZ 
DR.  D.  SIMONS 
DR.  G.  SARASCH 
DR.  L.  DUNCAN 
DR.  P.  BERNHARDT 
DR.  S.P.  GARY 

UNIVERSITY  OF  MARYLAND 
COLLECE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 


UNIVERSITY  Or  TEXAS 
AT  DALLAS 

CENTER  FOR  RESEARCH  SCIENCES 
P.O.  30X  538 
RICHARDSON,  TX  73030 
DR.  R.  HEELIS 
DR.  W.  HANSON 
dr.  j.?.  McClure 

UTAH  STATE  UNIVERSITY 
4TH  AND  3TH  STREETS 
LOGAN,  UTAH  34322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 
DR.  J.  ST. -MAURICE 

PHYSICAL  RESEARCH  LABORATORY 
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